Bayesian estimation of one- parameter qubit gates 



Berihu Teklu 

Dipartimento di Fisica, Universita di Milano, 1-20133 Milano, Italy 

Stefano Olivares 

CNISM, UdR Milano Universita, 1-20133 Milano, Italy 
Dipartimento di Fisica, Universita di Milano, 1-20133 Milano, Italy 



Matteo G. A. Paris 

Dipartimento di Fisica, Universita di Milano, 1-20133 Milano, Italy 

CNISM, UdR Milano Universita, 1-20133 Milano, Italy 

Institute for Scientific Interchange Foundation, 1-10133 Torino, Italy 



Abstract. We address estimation of one-parameter unitary gates for qubit 
systems and seek for optimal probes and measurements. Single- and two-qubit 
probes are analyzed in details focusing on precision and stability of the estimation 
procedure. Bayesian inference is employed and compared with the ultimate 
quantum limits to precision, taking into account the biased nature of Bayes 
estimator in the non asymptotic regime. Besides, through the evaluation of the 
asymptotic a posteriori distribution for the gate parameter and the comparison 
with the results of Monte Carlo simulated experiments, we show that asymptotic 
optimality of Bayes estimator is actually achieved after a limited number of 
runs. The robustness of the estimation procedure against fluctuations of the 
measurement settings is investigated and the use of entanglement to improve the 
overall stability of the estimation scheme is also analyzed in some details. 



1. Introduction 

Let us consider a system prepared in a known quantum state which enters an apparatus 
performing an operation on the state. The evolution imposed by the apparatus 
depends on the value of some parameters and the experimenter is interested in the 
estimation of those parameters. A natural strategy to obtain the parameter is to 
detect the state at the output and infer the value of the parameters from the global 
sample coming from a number of repeated measurements. The optimization of this 
strategy, i.e. the choice of the best probe, measurement and data processing, is 
the subject of quantum parameter estimation, which is a relevant subject in the 
quantum characterization of states and operations [H [2] • The operation on the state 
may be unitary or not [3l |4] and may depend on one or more unknown parameters, 
which, in turn, may correspond to quantities that are not directly observable. The 
parameters of interest may be the amplitude of the carrier signal, the position and 
orientation of an object, or the strength of an external fields. Communications, 
image analysis and precision metrology provide relevant examples. In the simplest 
scenario, a parameter estimation problem consists in the determination of the value 
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of the interaction parameter 6 for unitaries of the form U$ — {— i8G} where G is a 
Hamiltonian operator that generates the transformation. 

Generally speaking, quantum estimation is concerned with the problem of 
finding optimal ways to estimate quantum states and processes. In turn, it has 
recently attracted much interest in quantum information [5j [3 E] as a tool for 
characterization of signal and gates at the quantum level [101 HU H21 E2 HH H3 [TBI [17l 
118) . The canonical way to address estimation of states and operation is by quantum 
tomography, (see [19] for a review) i.e by measuring a complete set (a quorom [14]) of 
observables, which allows or the complete characterization. For single- and two-qubit 
systems this involves the measurement of Pauli matrices and has been realized for 
polarization qubits [201 HI] ■ Process Tomography, i. e the reconstruction of quantum 
operations [22J, [23l [24], is itself critical for verifying the actions of quantum logic 
gates 25 and characterizing decoherence processes [26]. On the other hand, there 
are many situations where the full tomography of signals and devices is not needed, 
either because the focus is on specific features of the transformation, or the dynamics 
is partially known. In this cases the relevant point is to find an optimal and stable way 
to achieve quantum characterization by parameter estimation. For this reason, in this 
paper we address estimation of one-parameter unitary gates for qubit systems, i.e. 
transformation of the form Ug = {—i8G} where G is a combination of Pauli operators 
and 8 is the parameter of interest. We consider the gate probed either by one-qubit 
and two-qubit states and compare the performances of standard measurements with 
the ultimate quantum limit to precision (accuracy) of estimation. As we will see, 
ultimate bounds are determined by the initial quantum state of the probe, the type 
of interaction and the readout measurements that is used to extract information from 
the probe. In particular, we are going to assess the performances of Bayes estimators, 
which themselves play a central role in many signal processing problems [27] . 

The precision (variance) of any unbiased estimator of a parameter 8 is limited 
by the Cramer-Rao bound (CR), given by the inverse of the Fisher information 
[28l [29l [30l [311 E21 [33] . Bayes estimators are known to be asymptotically unbiased 
and, in turn, to saturate CR asymptotically. For measurements that are related to 
the unknown 8 through a linear Gaussian model, the maximum likelihood estimate 
of 8 also achieves the CR. Furthermore, when 8 is estimated from independent 
identically distributed (iid) measurements, under suitable regularity assumptions on 
the probability density, the maximum likelihood estimator is asymptotically unbiased 
and achieves the CR |29[ [34] , On the other hand, being interested in realistic 
measurement schemes, here we consider estimation procedure based on a limited 
number of measurements. As a consequence, we have to take into account the biased 
nature of Bayes estimators. The variance of any estimator with a given bias is bounded 
by the biased CR [35l [36] , which is an extension of the CR taking into account the 
a priori distribution of the parameter of interest. In turn, it is a fundamental rule of 
estimation theory that the use of prior knowledge leads to a more accurate estimator. 

In this paper we address the estimation of the interaction parameter of unitary 
qubit transformations. We derive ultimate quantum limits to precision and assess 
performances of Bayesian estimators 137, ,38]. In particular, we focus our attention 
on measurement schemes as those in Fig. Q]and Fig- where a single-qubit gate is 
probed by single- and two-qubit probes, respectively. We evaluate the a posteriori 
distribution for the gate parameter, derive the ultimate bound on precision, and 
compare the asymptotic performances of Bayes estimator to that of Monte Carlo 
simulated experiments, thus showing that asymptotic optimality is achieved after a 
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limited number of runs. The paper is structured as follows. In Section [5] we introduce 
notation and derive the a posteriori distribution, also discussing the Bayesian version 
of the CRB. In Section [3] we discuss limits to precision in estimating unitary gates for 
qubit systems. A comparison between single- and two-qubit entangled probes shows 
that entanglement improves the overall stability of the estimation procedure. We also 
compare the asymptotic a posteriori distribution for the gate parameter to the results 
of the Monte Carlo simulated experiments. Section [4] closes the paper with some 
concluding remarks. 

2. Parameter estimation of one-parameter qubit gates 

A generic unitary transformation acting on a qubit state can be written as U(6) = 
exp{— i;0 ■ <x}, where a = {<Ti, (T2, er 3 } is the vector of Pauli matrices and is a vector 
describing the transformation. In the following we will assume that depends on a 
single parameter and refer to it as to 9. The simplest scheme for the estimation of 9 
[1] consists in choosing a probe pure qubit state |^o) undergoing the transformation 
and a measurement onto the evolved state \ipg) = U(0)\ipo). Here we assume that 
the measurement can be represented by the two projectors Hq and Hi, so that 
the conditional probabilities to obtain the outcomes "0" or "1" (i.e. given 0) are 
P(j\6) = (ipo\ILj\ipo) , j = 0, 1. After M measurements on equally prepared qubits we 
have the the sample X = {x\, . . . , xm} of outcomes, where the x^s can take the values 
"0" and "1". This leads us to define the following sample probability or likelihood 
function: 

M 

P(X\6) = l[P(x k \d). (1) 

k=l 

Our estimation problem is that of inferring the value of 9 once the sample of outcomes 
is assigned by the measurement; in other words, we are interested in the conditional (o 
posteriori) probability P(9\X) of 9 given the sample X. This can be easily obtained 
by the Bayes theorem, which states that [39j [40] P{9\X)P(9) = P(X\9)P(X) where 
P{9) is the prior probability and P(X) is the overall (unconditional) probability of 
the observed sample. Hence, the a posteriori distribution may be written as 

1 M M 

P{9\X) = -\{P{x k \9), N= d9\{P{x k \9), (2) 

k=l ^° k=l 

where N is the normalization and Vl is the set of possible values for 9. From ([2]) we 
may evaluate the expected value of 9 and the variance of the distribution 

9= fd99P{9\X), Var[0]= (d9{9-9) 2 P{0\X) . (3) 
in J q 

The mean (expected) value 9 of the a posteriori distribution is our Bayesian estimator. 

For a large number of measurements, M ^> 1, and assuming that the true value 
of the parameter is 9* , the number of times a factor P(x\9), with x = 0, 1, appears 
in the product © is approximately given by P(x\9*)M . The asymptotic a posteriori 
distribution for the parameter 9, conditioned on the true value 9* is thus given by 

mum 

n pm p{sinM - (4) 
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Since E/,=o,i p ( h \ d ) = 1 we have d e P M (0\0*) \ g , = and d%P M (9\0*) \ g . < 0, i.e the 
distribution Pm(9\9*) has the desirable property of showing a maximum at the true 
value of the parameter, i.e. Bayesian estimator is asymptotically unbiased. 

According to the Laplace-Bernstein-von Mises theorem [33] the a posteriori 
distribution of Eq. Q may be asymptotically approximated by a Gaussian with 
variance given by a 2 = [MG(e*)Y 1 where we have introduced the Fisher information 



_ v 1 f dP(s\9) 



(5) 



The asymptotic a posteriori distribution is thus completely characterized by its 
variance, or equivalently by the Fisher information which itself gives a lower bound 
to the variance of any unbiased estimator 9(X) via the Cramer-Rao inequality 
0[3J[3J[33]: 

Any estimator saturating the inequality © is referred to as an efficient estimator. The 
relation a 1 = [MG(8*)] 1 thus says that Bayesian estimator is asymptotically efficient, 
whereas this conclusion does not hold for finite M. For finite M Bayes estimator is 
biased and we have to generalize Eq. ((6]) to take into account the bias. To this aim 
one considers the conditional expectation of the error B(9) = J dx [9(x) — 9]P(x, 9), 
where P(x,6) — P(x\9)P(0) is the joint probability of the data and the parameter. 
Of course, for unbiased estimators B(9) = 0. Starting from the definition of B(9) one 
derives the so-called van Trees inequality 35j for the mean squared error 



Var[0(x) -8]= J d8 P{9) J dx [8{x) - 0\*P{x\9) > , (7) 

where we introduced the generalized Fisher information Hm{&) — F(9) + MG(ff), G 
being the Fisher information of Eq. |(SJ), M the number of repeated measurements, 
and F the Fisher information of the prior, i.e. F{9) = f d6 [d e log P{8)} 2 P{9). Eq. © 
takes into account the information due to the prior and thus gives a lower bound than 
the CR one, which is anyway achieved for M ^S> 1. On the other hand, one may show 
that [36]: 

H M {8)= JdXP(X) Jd9[d e \ogP{8\X)] 2 P{8\X) 

M ^ F M {8\8*) ee Jd9P M {8\8*) [d e log P M (8\8*)] 2 , (8) 

where P(9\X) is the Bayesian probability distribution of Eq. ([2]) and Pm(9\9*) its 
asymptotic expression of Eq. ((4]). Thus, Eq. (© represents the Bayesian counterpart 
of the CR. 



3. Bayesian estimation of one-parameter qubit gates 

In this section, without loss of generality, we address the case in which = (0,0,9), 
9 6 [0, 7r], i.e., 1/3(9) = exp(— ^9cr 3 ). We first consider the gate probed by a single- 
qubit state and jointly optimize the probe and the measurement and then address the 
use of entanglement, showing that it may be useful to improve the overall stability of 
the estimation procedure. 
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3.1. Estimation via sinqle-qubit probes 

In Fig.Q]we schematically depict the estimation procedure: a pure state \tpo) undergoes 
the unitary transformation 1/3(6) and, then, is measured by means of a projective two- 
outcome device. Our aim is that of optimizing the estimation of 6 by a suitable choice 
of both the probe state and the projective measurement. 





U 3 (9) 




► 


>■ 



Figure 1. Estimation of qubit gates. The state {po} is prepared, then enters the 
gate described by the unitary U%(6) with unknown 6, and finally is detected by a 
detector described by a two-value POVM. 

Upon writing the pure qubit state in the standard way 

^o)=cos||0)+e^sin||l), (9) 

where |0) and |1) are the two eigenvectors of <r 3 , the parameters a <G [0, 7r] and 
4> € [0, 2tt] uniquely determine |V>o) and the evolution under the unitary 1/^(9) is 
straightforward. Next, we perform the measurement described by the two projectors 

n (M = ni09 > w) = i-noC9,w) ) (10) 

where \V(j3,u)) = cos||0) + e iuJ sin ^ 1 1} . The probabilities of the two outcomes are 
given by P(0\6) = P (a, (3,0^,6) = \ (0\U e \ft>o) | 2 ie 

P(0\9) = — [1 + cos a cos (3 — cos(0 — oj + 9) sin a sin (3] 
P(l\6) = l-P(0\9). (11) 
The Bayesian a posteriori distribution of Eq. @ may be written as 

P(0\M) = j- P{0\9) rno P(l\9) mi , (12) 

where rnj is the number of measurements with outcomes j = 0, 1, mo + mi = M , in 
the observed sample. For a large number of measurements M 3> 1, we can evaluate 
Pm{9\9*) using Eq. (|4]) and, in turn, the expectation 9 and the variance Var[0]. Upon 
expanding on a and j3 up to second order one sees that Var[0] achieves its minimum 
for the choice a = (3 = ir/2, independently on the value of </> and ui. This is confirmed 
by the evaluation of the corresponding Fisher information G(6), 



sin a sin /3sin 
1 — cos 2 a cos 2 (3 — cos 2 (9 + <ft — uS) sin^ a sin^ f3 



Bin usiii /Join ^Ti?-^ 
G W = i ~ 9 a .._9/» , I . A Z 2 ■ 2 o ' 



which achieves its maximum G = 1 for a = f3 = 7r/2, independently on the value of > 
and W. Using this results we have: 



M ( cos 2 — log cos 2 — h sin 2 — log sin 2 - 
1 2 2 2 2 



(14) 
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and the variance saturates the van Trees inequality, thus confirming that Bayes 
estimator is asymptotically efficient. In Fig. [5] we report the ratio 8/8* and variance 
multiplied by Hm [see Eq. 0] as a function of the number of measurements. As it 
apparent from the plots all the curves approach one when the number of measurements 
increases. In the asymptotic region Hm — MG(8) ~ M. 

8/8* 
1.003K 

1.002 
1.001: 



1 .000 [ . = 

300 500 700 

Figure 2. LogLinear plot of 6/6* (left) and linear plot of ~Vai[9]H m (right) as 
a function of the number of measurements M for the estimation of the unitary 
Uz{9) and for different values of the true parameter. For both plots the dotted- 
dashed line is for 9* = 0.8, the dashed line is for 9* = 1.2 and the solid line is for 
9* = 1.8. 




Notice that the results here reported for U^{8) are actually valid for any other 
unitary gate. This can be seen as follows: any unitary U n {8) — exp (—icr n 8) describing 
rotations around the arbitrary axis n may be written as U n — 011^,(8)0^ where O is 
the rotation corresponding to the mapping 3 — > n. As a consequence, the optimal 
measurement corresponds to the projectors IIq = 0+1^(0, w))(^'(0, u))\0, 11^ = 1 — IIq 
and the optimal probe is given by O^|0). 

A question arises on whether the results reported above may be easily 
implemented in practice. This concerns the stability of the measurement rather than 
its precision. The point is the following: Suppose that for some reasons the values of 
parameters a (probe) f3 (measurement) slightly deviate from the optimal settings. To 
which extent the overall performances of the procedure are degraded? A convenient 
way to address this issue is to make a perturbation analysis upon expanding the Fisher 
information G{8) around the optimal settings a = (3 = 7r/2 up to second order 

I cos 

G{8) ~ 1 - [(« - tt/2) 2 + 08 - n/2f] + 2^-(a - tt/2)(/3 - vr/2) 

sin 8 sin 8 

^ X l-i(a-/3) 2 . (15) 



Eq. (|15[) shows the quadratic decrease of G out of the optimal setting and, 
especially for small values of the gate parameter 8, the dramatic effect of a mismatch 
between the values of a and j3. The latter is well illustrated in Fig. [3] where G as a 
function of a and (3 is shown for different values of 8. As it is apparent from the plots a 
mismatch \a — f3\ ~ 8 of the order of the parameter to be estimated is enough to make 
the whole procedure ineffective. Fortunately, as we will see in Section |3~31 the stability 
issue may be overcome by using entangled probes and the optimal performances still 
achieved by a two-qubit probe configuration. 
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Figure 3. Contour plot of the Fisher information G(9) as a function of the 
probe and measurement parameters a and /3 for three different values of the gate 
parameter 9* = 0.05,0.1, 1. Darker areas corresponds to higher values. 



3.2. Monte Carlo simulated experiments 



Before addressing stability of the measurement let us compare the asymptotic a 
posteriori distribution for the gate parameter with the results of Monte Carlo simulated 
experiments. This is in order to locate the asymptotic region and make quantitative 
statements on the achievability of the ultimate bounds to precision. We have simulated 
M repeated measurements using the optimal probe/measurement settings and inserted 
the resulting values of m,j 1 j — 0,1 in Eq. (|12|) to obtain the a posteriori distribution 
for the gate parameter. In Fig. [4] we plot the rescaled variance of this distribution 
together with the variance of the asymptotic a posteriori distribution of Eq. (fnj. 
Remarkably the asymptotic region is achieved after a limited number of runs, thus 
proving that Bayesian approach may useful in practical applications. In the inset we 
report the full distribution, both the experimental a posteriori and the asymptotic 
one, for M = 20 and M = 500. Notice that: i) for M = 500 the a posteriori 
experimental distribution is already indistinguishable from the asymptotic one and 
ii) the asymptotic is already unbiased, i.e efficiency has been achieved for a limited 
number of runs. 




200 



400 



600 



Figure 4. Rescaled variance (variance multiplied by the generalized Fisher 
information) of the a posteriori distribution as a function of the number of 
measurements. Black line is for the asymptotic distribution whereas gray squares 
are for the experimental one. Inset: asymptotic (dashed line) and experimental 
(solid line) Bayesian a posteriori distributions for 9* = 0.6 as obtained using 
optimal single-qubit probe and M = 20 (left) or M = 500 (right) repeated 
measurements. 
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3.3. Estimation via two-qubit entangled probes 

An alternative scheme for gate estimation may be designed using entangled states 
as depicted in Fig. [5] In fact, the use of entanglement may improve estimation 
[44l l45l l46l |47] . In this section we investigate whether this is the case for the present 
estimation problem. Basically, the use of entanglement increases the dimension of 
Hilbert space and thus the number of possible outcomes of a measurement performed 
on the perturbed signal. The corresponding Fisher information does not increase 
but the maximum value is achieved for a large class of probe signals. The Bayesian 
estimator is able to exploit this fact to increase the overall stability of the estimation 
procedure of the gate parameter. 




■ H^ (e) I — 





k= 0,1,2,3 



Figure 5. Schematic diagram of parameter estimation by entangled qubit probe: 
{IV'o))} is prepared and is subjected to the unitary transformation Ug on one 
qubit. Finally, a projective measurement {n a } is performed. 

In order to see this behavior in practice let us consider the estimation of the 
parameter of the gate U^{9) using a generic pure state of the form 

l^>> = ^X>k*>> ( 16 ) 

where we used the matrix notation for states \A)) = J^ij Ay = A® As 
a measurement we consider the "Bell" measurement made of the four projectors 
life = !|(7fe))((o-fe|, fc = 0,1,2,3, over a set of maximally entangled states. After 
the evolution under the unitary U^{6) the four possible outcomes of the measurements 
occur with the probabilities 

9 6 {06 

P Q (9)= cl cos 2 - + c 2 sin 2 - , P x {6) = I ci cos - - c 2 sin - 

P 3 (6)= cl sin 2 9 - + c 2 3 cos 2 t , P 2 (6) = U sin °- + c 2 cos ~ ) ( I 7 ) 

The corresponding Fisher information is given by 

(cl-cj)(4~cj)sm 2 6 



G(9) =ct + 4 + 



( c o + c §)-( c u- c !) cos2 



and achieves its maximum [G = 1, \/9) for any state with cq = or C3 = 0. Therefore, 
having fixed the Bell measurement at the output, a possible deviation in the probe 
preparations is not degrading the performances of the estimation procedure. 

4. Conclusions 



In this paper we have analyzed estimation of one-parameter unitary gates for qubit 
systems. We have addressed Bayesian estimation procedures and compared their 



Bayesian estimation of one-parameter qubit gates 



performances with the ultimate quantum limits to precision. Bayes estimator is 
known to be asymptotically unbiased, but for practical implementation is of interest to 
evaluate quantitatively how many measurements are needed to achieve the asymptotic 
region. To this aim, after the evaluation of the asymptotic a posteriori distribution 
for the gate parameter, we have compared it to the distribution obtained by Monte 
Carlo simulated experiments and full Bayesian analysis and shown that asymptotic 
optimality of Bayes estimator is achieved after a limited number of runs. We have 
also addressed the issue of stability, i. e the robustness of the optimal settings against 
fluctuations of the probe and measurement parameters. It has been shown that the 
use of entanglement is useful to improve stability. More explicitly, we have shown that, 
although the Fisher information does not increase, its maximum value is achieved for 
a large class of probe signals, thus making the procedure more robust and increasing 
the overall stability of the estimation procedure. 
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